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INTRODUCTION 

The r\ decay channels with lepton pairs are closely related to the channels with real 
photons. The direct consequence of Quantum Electrodynamics is that the process with 
a real photon should be accompanied by a process where a virtual photon convert 
internally into a lepton pair. This fact was first pointed out by Dalitz in 1951 yj] and 
therefore the decays with one or two lepton-antilepton pairs are known as single or 
double Dalitz decays. Dalitz decays can be related to the corresponding radiative decays 
using Quantum Electrodynamics and introducing a function of the four momentum 
transfer squared of the virtual photon (q 2 ): F(q\,q 2 ) - the Form Factor. In case of 
Dalitz decays q 2 is equal to the invariant mass squared of a lepton pair and q 2 > 4-m 2 . The 
Form Factor describes the structure of the transition region and it is also used to describe 
processes where q 2 < (space-like virtual photons). The orientation of the plane defined 
by the lepton momenta is connected with the polarization of the virtual photon. 
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FIGURE 1. Diagrams for a) single and b) double Dalitz decays of the pseudoscalar meson rj. Diagram 
c) is for the dominating conventional mechanism for rj meson decay into lepton-antilepton pair. 

The diagrams corresponding to the conversion decays of r\ are shown in fig. dkb. 
The diagram for the very rare decay into a lepton-antilepton pair P — > presented in 
fig- QJ, is forbidden to proceed via a single photon intermediate state. The dominant 



mechanism within the Standard Model is the second order electromagnetic process 
with two virtual photons P —> y*Y*- m addition the decay is suppressed by helicity 
conservation. Therefore it is a potentially important channel to look for effects of physics 
beyond the Standard Model. The imaginary part of the decay amplitude can be exactly 
related to the decay 77 — > 77. The measured value of r(rj — »• 77) gives a lower limit (the 
unitarity limit) of BR(T] — > e+e~)>1.8x 10 -9 M,\M which is a much lower value than 
for other decays of n° and 77 into lepton-antilepton pairs. This makes the 77 — > e + e~ 
decay sensitive to contributions from non-conventional effects that would lead to a 
higher branching ratio. The current experimental upper limit of 7.7 x 10~ 5 (90% CL) 
obtained by the CLEO II collaboration [|4J] is four orders of magnitude higher. The real 
part of the amplitude of the 77 — > e + e decay can be estimated using the measured value 
of BR(rj — > S 0, ED and gives a prediction for the branching ratio of about 

6 x 10 . The estimated BR is based on the assumption that the ratio between Im and 
Re part of the decay amplitude is the same for e + e and jU + /i . 

Dalitz decays of the 77 allow studies of the structure of the 77 — > 7*7 and 77 — > 7*7* 
vertex. An interesting question is how the Form Factor F(q\,q^) behaves as a function 
of the invariant mass squared of the pairs q 2 . Different forms of this dependence have 
been proposed, the most common is predicted by the vector meson dominance model. 
Experimental information for the 77 Dalitz decays are scarce. The branching ratio for 
77 — > e + e~y decay is known with rather large error BR=(6.0±0.8)xl0~ 3 [8]. The 
rate of the double lepton pair decays of the 77 has not been determined experimentally, 
but different theoretical predictions can be found in IsLQll [HI [12] • The measured and 
calculated branching ratios for various lepton decay channels are summarised in tableQ] 



TABLE 1. The measured and calculated Branching Ratios for different 
77 decay channels with lepton pair(s). 



Decay mode 


Branching Ratio exp. 


Branching Ratio theor. 


1. 77 — > e + e~y 


(6.0 ±0.8) x 10- 3 


(6.37-6.57) x 10~ 3 


2. rj^n+n-y 


(3.1 ±0.4) x 10~ 4 


(2.1-3.05) x 10~ 4 


3. 77 — y e + e~e + e~ 


<6.9x 10~ 5 


(2.52 - 2.64) x 10~ 5 


4. 7] — y e + e~H + H~ 




(1.57-2.21) x 10~ 7 


5. 77 -> n+n'ii+ii- 




2.4 x 10~ 9 


6. 77 > e+e~ 


<7.7x 10~ 5 


1.8 x 10~ 9 1 


7. T]^fl+fl- 


(5.8 ±0.8) x 10~ 6 


4.4 x 10- 6 ~ 


8. 77 — ► n + n~e + e~ 


(4.3 ±1.3 ±0.4) x 10~ 4 


3.1 x 10~ 4 


9. 77 > n + 7t~jj, + ii~ 




7.5 x 10~ 9 


* unitary bound 



THE EXPERIMENT 

The experiment was performed at the CELSIUS storage ring in Uppsala, using the 
WAS A detector setup (fig. [2} lfl3ll . Beam protons with a kinetic energy of 893 MeV 
interacted with frozen droplets of deuterium [14]. The 77 mesons were produced in the 
reaction pd ^ 3 He 77 close to the 77 production threshold. The detection of 3 He ions in 




a zero-degree spectrometer (tagging detector) provided a clean r\ trigger independent of 
decay channel. The 3 He ions were identified and their energy was measured which al- 
lowed to select cleanly the pd — > 3 He r\ reaction with a background of about 1%. The tag- 
ging detector provided a few triggers per second (at a luminosity of 5 x 10 cm s ), 
yielding on average one recorded r\ event per second. During two weeks of the experi- 
ment (distributed over half a year period) nearly 3 x 10 5 r\ events were collected. 

The charged 77 decay products are tracked using a cylindrical mini drift chamber 
(MDC), consisting of 17 layers of thin mylar tubes and built around a beryllium beam 
pipe with 1.2 mm wall thickness. The MDC is placed inside of a superconducting 
solenoid which provides a magnetic field of IT The MDC is surrounded by a barrel of 
plastic scintillators. An electromagnetic calorimeter consisting of 1012 CsI(Na) crystals 
measures the energies of photons and impact points. Correlation between momentum, 
measured from the track curvature in the MDC and energy deposit in the electromagnetic 
calorimeter permits to distinguish electrons from charged pions if the transverse particle 
momenta are larger than about 20 MeV/c. 

In the offline analysis, events with at least two charged particle tracks reconstructed 
in the MDC were required to originate from a point close to the beam target intersection 
region. Hit clusters in the calorimeter, without associated tracks in the MDC and with 
energy deposit larger than 20 MeV were assumed to originate from photons. Only 
events containing decay particle candidates with balanced electric charge were accepted 
for further analysis. The results on the rj — > Tl + ll~e + e~ decay channel were already 
presented in ref. [fl5|l . Events with a pair of charged decay products with opposite electric 
charges can be attributed either to the decay channels with two charged leptons or to 
more frequent channels with two charged pions (ri — > n + n~y and 77 — > 7C + n~n°). 

The following variables were used in the further data analysis: 

1. The invariant mass of the pair of charged particles M ee . The electron masses were 
assumed in the calculations. A clear peak at the smallest possible mass value is 
expected for e + e Dalitz pairs and from conversion of real photons in the detector 
material. 



2. The invariant mass for all decay products. It was required to be in the region of the 
r\ mass. 

3. The missing mass of the total decay product system (MMrj)- It should be, within 
errors, equal to the mass of the 3 He nucleus. 

4. The ratio between the momentum measured in the drift chamber and the energy 
of the shower in the electromagnetic calorimeter associated with the charged track 
Rp/E - It permits to distinguish between electrons and charged pions for the particles 
reaching the electromagnetic calorimeter. 

NORMALIZATION 

In order to normalize the branching ratios of the r\ meson decays involving an electron- 
positron pair a monitoring process is needed to check the reconstruction efficiency 
for the electrons. Dalitz decay of one of the three n° from the r\ — > 11®%®%® decays 
provides an abundant data set of events with five photons and an electron-positron 
pair - r\ — > 7C°n°np. The Dalitz decay of the n° meson has been studied in detail 
both theoretically and experimentally. According to the Monte Carlo (MC) simulations 
already the class of events with more than three neutral clusters is dominated by this 
decay channel. The invariant mass of the two charged particles for such events is 
presented in fig. [3^. The peak at low masses can be attributed to the e + e~ pairs from 
the neutral pion decay. The expected distribution from MC is shown in the same figure 
with dashed line. The second maximum is due to the decay channels with n + iz~ pair, 
mainly r\ — ► (the corresponding histogram in the figure is normalized to the 

second enhancement). Such interpretation is confirmed by the reconstructed invariant 
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FIGURE 3. Event sample with two tracks in MDC and at least four neutral clusters - candidates for 
77 — > K°n°np decay. Distributions of: a) M ee - experimental data (points), MC for rj — * n Q n np decay 
(dashed line), MC for decays involving a misidentified n + n~ pair (dotted histogram), b) M ee y - all data 
(solid line) and events with M ee < 0.1 GeV/c 2 (dashed line). 

mass of the e + e~y system (M ee y) where the photon leading to the mass value closest to 



the n° mass is selected (fig. [3b). The distribution leads to a peak at n° mass when the 
M ee < 0. 1 GeV/c 2 condition is applied. 



RESULTS 



The candidates for r\ — > e + e~y decay are searched in the class of events with less than 
four neutral clusters. One of the clusters should have energy deposit greater than 200 
MeV. More than one cluster is allowed in order not to limit the acceptance for the 
r\ — > e + e~y decay, where due to electron or photon interaction in the calorimeter an 
additional hit cluster can be created. In fig. |4]the experimental M ee distribution is shown 
together with MC distributions for the r\ — > e + e~y decay and a background from rj 
decays with a pair misinterpreted as electrons. The plot is before identification 

using conditions on the R p /E variable. 
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FIGURE 4. M ee distribution for events with less than four neutral clusters - the initial selection for 
Tj — » e^e~y decay. Histogram - experimental data, filled area - MC simulation for rj — > e+e~y decay 
and background from decays involving a misidentified K^K~ pair. 
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FIGURE 5. Selected data sample of the T] — > e^e~y decay candidates: relative angle between momen- 
tum of the electron-positron pair and the photon momentum vs M ee y invariant mass, (a) data, (b) MC 
simulations for r\ —* e+e~y. 
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In the final selection of the 77 — > e + e~y decay candidates it was required additionally 
that at least one electron is identified using R p /E rat i° an d M ee < 0. 125 GeV/c 2 . In fig. [5] 
the relative angle between the momentum of the electron-positron pair and the photon 
momentum is plotted with respect to the invariant mass M ee y after the selection of the 
r\ — > e + e~y candidates. Fig. [6] shows the experimental M ee y distribution (points) and 
the result of the MC simulation of the signal and the background (solid line) after the 
selection cuts. The origin of the discrepancy between MC and the experimental data in 
the figure around 0.42 GeV/c 2 is under study. It can be due to the interaction of the beam 
with rest gas in the beam tube. 
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FIGURE 6. Selection of the rj — * e + e y decay candidates. The invariant mass M ee y distribution: 
experimental data (points); sum of the MC simulation of the 77 — * e + e~y and the background channels 
(line). 

In a search for the 77 — ► e + e~e + e~ decay, events with exactly four tracks in the MDC 
and equal number of positive and negative charged particles were selected. According 
to the simulations 11% of the reconstructed 77 — > e + e~e + e~ events fulfill the following 
criteria: 

1 . the relative angle between electron and positron in both pairs is less than 40° 

2. the relative angle between the momenta of the two e + e~ pairs is in the interval 
(110°-170°) 

3. the 77 meson emission angle is smaller than 45° 

4. the missing transverse momentum is less than 0.3 GeV/c 

5. the missing mass MMfj is in the interval 2.65-2.90 GeV/c 2 . 

In the data only two events passed all selection cuts. A schematic view of the event 
display for one of the two candidates is shown in fig.[7£Left). The estimated background 
is 1.0 ± 0.3 events mainly from 77 — > e + e~y with a y conversion into e + e~ pair in the 
detector material. 

The v\ — > e + e~ decay has a distinctive signature in the pd ^ 3 He r\ reaction close 
to threshold: the emitted electron and positron have large energies (E > 150 MeV), are 
in plane with the beam and large opening angle (about 134°). In fig. [TtRight) the e + e~ 
invariant mass is plotted as a function of the opening angle of the electrons for the 
data sample and for the MC simulation of the 77 — > e + e~ decay. The selection included 
the electron identification condition: 0.5 < R p /E < 1-65. No events are observed in 
77 — > e + e~ signal region determined by MC simulations (corresponding to the overall 
reconstruction efficiency of 41.2%). 
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FIGURE 7. (Left) Event display for an rj — > e + e~e + e~ candidate event. The shaded area in the 
outermost ring represents the projection of the hit calorimeter crystals (the size of the crystals and 
the radial position of the front faces is not to the scale). (Right) Relative angle between electron and 
positron tracks vs M ee for tj — > e + e~ event sample selection: data - scatter plot; MC simulation of the 
tj — > e + e~ decay - shaded area. 



SUMMARY 

A data sample of about five hundred r\ — > e + e~y decays events were collected. This 
number is consistent with the normalization based on r\ — > tz q tz q tZq decay events. Two 
candidate events for the double Dalitz decay mode r\ — > e + e~e + e~ were observed with 
an estimated background of 1.0 ±0.3 events. It is the first observation of this r\ decay 
channel. No candidates have been found for the r\ — > e + e~ decay in a sample of about 
230k clean r\ events. 
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